Experimental investigation of the dynamics of entanglement: Sudden death, 
complementarity, and continuous monitoring of the environment 
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We report on an experimental investigation of the dynamics of entanglement between a single 
qubit and its environment, as well as for pairs of qubits interacting independently with individual 
environments, using photons obtained from parametric down-conversion. The qubits are encoded in 
the polarizations of single photons, while the interaction with the environment is implemented by 
coupling the polarization of each photon with its momentum. A convenient Sagnac interferometer 
allows for the implementation of several decoherence channels and for the continuous monitoring of 
the environment. For an initially-entangled photon pair, one observes the vanishing of entanglement 
before coherence disappears. For a single qubit interacting with an environment, the dynamics of 
complementarity relations connecting single-qubit properties and its entanglement with the environ- 
ment is experimentally determined. The evolution of a single qubit under continuous monitoring of 
the environment is investigated, demonstrating that a qubit may decay even when the environment 
is found in the unexcited state. This implies that entanglement can be increased by local continuous 
monitoring, which is equivalent to entanglement distillation. We also present a detailed analysis 
of the transfer of entanglement from the two-qubit system to the two corresponding environments, 
between which entanglement may suddenly appear, and show instances for which no entanglement 
is created between dephasing environments, nor between each of them and the corresponding qubit: 
the initial two-qubit entanglement gets transformed into legitimate multiqubit entanglement of the 
Greenberger-Horne-Zeilinger (GHZ) type. 

PACS numbers: 03.65.Yz; 03.67.Bg; 03.67.Mn; 42.50.Ex 



I. INTRODUCTION 

Entanglement plays a central role in quantum me- 
chanics. The subtleties of this phenomenon were first 
brought to light by the seminal paper of Einstein, Podol- 
ski, and Rosen JE, published in 1935, and by those of 
Schrodinger 0, [J published in 1935 and 1936. It took 
however approximately thirty years for its essential dis- 
tinction from classical physics to be unmasked by John 
Bell [4|, and another thirty years for the discovery that 
entanglement is a powerful resource for quantum com- 
munication 0, |{| @, H, H, Qlj ■ It was also found in the 
90's to play an important role in quantum computation 
algorithms [ilj . Furthermore, it plays a key role in the 
behavior of macroscopic quantities like the magnetic sus- 
ceptibility at low temperatures 12 1. 

Yet the dynamics of entangled systems under the un- 
avoidable effect of the environment is still a largely un- 
known subject, in spite of its fundamental importance 
in the understanding of the quantum-classical transition, 
and its practical relevance for the realization of quantum 
computers. 

The absence of coherent superpositions of classically 



distinct states of a macroscopic object is analyzed by 
decoherence theory [l3|, [l4| , which shows that the emer- 
gence of the classical world is intimately related to the 
extremely small decoherence time scale for macroscopic 
objects. Within a very short time, which decreases with 
the size of the system, an initial coherent superposition 
of two classically distinct states gets transformed into a 
mixture, due to the entanglement of the system with the 
environment. The decay dynamics is ruled, within a very 
good approximation, by an exponential law. 

Detailed consideration of the dynamics of entangled 
states requires defining proper measures of this quan- 
tity. For pure states, one can use the Von Neumann 
entropy [ll| associated to each part, or alternatively the 
corres pon ding purity, defined by the so-called linear en- 
tropy [l5l . UU . The ideais that the more entangled some 
partition of a multiqubit state is, the more unknown is 
the state of each part. 

However, systems undergoing decoherence do not re- 
main pure. A mixed state of N parties is separable if it 
can be written as a convex sum of pro ducts of density 
matrices corresponding to each part [171 ] : 
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where the index /i refers to the /x-th realization of the 
state and = I; with > 0. 
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Entanglement measures for mixed states have been de- 
fined for systems with dimension up to six [IE [IE H(| , 
but for larger dimensions this problem has not yet been 
solved. For two-qubit systems, Wootters [I|| introduced 
the concurrence as a measure of entanglement. 

It was shown by Peres [l9| that, if the partial trans- 
pose of the density matrix of a multipartite system with 
respect to one of its parts has negative eigenvalues, then 
the state is necessarily entangled. Thus, a non-negative 
partial transpose is a necessary condition for a state to 
be separable. However, this condition is also sufficient 
only for 2 x 2 or 2 x 3 systems, as shown in Ref. [20j . 

The negativity, defined as the magnitude of the sum 
of negative eigenvalues of the partially transposed ma- 
trix, can thus be used, in these cases, as a measure of 
entanglement. For higher dimensions, this does not work 
anymore: the negativity is then an indicator of distill- 
able entanglement. That is, if the negativity is different 
from zero then it is possible, through local operations 
and classical communication, to obtain from n copies of 
the state a number m (m < n) of maximally entangled 
states. This process, called distillation, does not work if 
the partially-transposed density matrix is non-negative: 
any entanglement still present cannot be distilled - it is 
then called bound entanglement [2~i| . 

These measures allow one to study the dynamics of 
initially entangled states under the influence of the envi- 
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Hi], HE l36j|. The outcome of these investigations is that 
the dynamics of entanglement can be quite different from 
the dynamics of decohcrence: the first is not ruled by an 
exponential decay law, , as the latter, and entanglement 
can disappear at finite times, even when system coher- 
ences decay asymptotically in time. This phenomenon, 
known as entanglement sudden death [32| . is a peculiar 
feature of global dynamics. 

In this article we present an all-optical device to study 
the interaction of simple systems (one or two qubits) with 
various kinds of environments, in a highly controllable 
fashion. The setup is extremely versatile, allowing to im- 
plement many different types of open system dynamics. 
A partial account of our experimental results was given 
in Ref. [33|. Here we show how this set up can be used 
not only to demonstrate the subtle dynamics of entangle- 
ment, but also the behavior of a continuously-monitored 
system, as well as the dynamics of complementarity rela- 
tions [13, HE HE HO] between local and global properties 
for a two-qubit entangled system. These complementar- 
ity relations quantify the notion that, for pure entangled 
states, coherences and populations of each party become 
uncertain: the more unknown they are, the more entan- 
gled is the state. 

Sections |TT] and IIIII contain the theoretical framework, 
in a form which is particularly suitable to the experimen- 
tal investigation of the dynamics of entanglement under 
different kinds of environment. 

Section |TT] deals with open-system dynamics and Kraus 
operators, while Section [Ml discusses quantum channels, 



the dynamics of complementarity for each of these chan- 
nels, and the transfer of entanglement from the two-qubit 
system to the two corresponding environments. A pecu- 
liar feature of dephasing processes is emphasized, for a 
family of initial states: when the two-qubit entanglement 
disappears, no bipartite entanglement is left in the sys- 
tem. The state of the two-qubit system plus correspond- 
ing environments becomes a state of GHZ (Greenberger- 
Horne-Zeilinger) type (4lj , with only genuine multiparti- 
cle entanglement. 

The experimental setup is introduced in section IIV1 
along with several examples of environments that we are 
able to implement. In section [V] we present the exper- 
imental results for the behavior of a single qubit, with 
and without continuous monitoring of the environment, 
including a detailed study of the dynamical behavior of 
complementarity relations between local (single-party) 
and global properties of the system qubit+environment. 
We show that our results on the continuously-monitored 
system are intimately related to the distillation of entan- 
glement. 

In Section lVTl we discuss the experimental investigation 
of the evolution of entanglement for two typical noise 
channels - amplitude damping and dephasing - including 
the first observation of the phenomenon of entanglement 
sudden death, which we had previously reported in [33| . 
Our conclusions are summarized in section IYIII 



II. OPEN SYSTEM DYNAMICS AND KRAUS 
OPERATORS 

A system (S) interacting with an environment (E) is 
described by the following Hamiltonian: 

H = H S (8)1 + 1® He + XVse, (2) 

where H$ and He are the system and environment 
Hamiltonians respectively, and Vse is the coupling term 
between them with coupling constant A (in the weak cou- 
pling limit, Ac 1). The system and the environment get 
entangled due to the interaction Vse ~ an initially pure 
state of S evolves to a mixed state. 

In quantum optics, the traditional way of dealing with 
open systems weakly coupled to environments with large 
number of degrees of freedom is through master equa- 
tions (421 . |43[ . In this approach, the equation of motion 
for the state ps of the system, given by: 

2 

PS = --Tr E [H,p SE ), (3) 

where pse is S + E density matrix, is approximated to 
first order of perturbation theory, with additional as- 
sumptions of Markov dynamics and initially uncorrelated 
systems. The previous expression can be then written as 
a sum of a unitary contribution plus a non-unitary term, 
which depends only on operators acting on the system S, 
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and is given by the following expression: 



•NU 

Ps 



CiCkps - 2C k psC\ 



(4) 



where the upper index NU stands for non-unitary, and Ck 
are the so-called Lindblad operators. See References |44l . 
|45| for a comprehensive treatment. 

The experimental investigation of open system dynam- 
ics can be greatly simplified by adopting an alternative 
formalism, based on the Kraus representation (46[. We 
summarize in the following the main ingredients of this 
approach. 



A. Kraus operators 

As suggested in Eq. ([3]), the evolution of a system cou- 
pled to an environment can always be expressed as a uni- 
tary dynamics on a higher dimensional system - Fig. [1] 
depicts this approach. 

Starting with uncorrelated systems, the total evolution 
can be written as: 



U S e(ps ® \0) E (0\)Ul 



SE i 



(5) 



where Use is the S + E evolution operator and |0) B , 
without loss of generality, represents the initial state of 
the environment. If we wish to focus only on the evolu- 
tion of system S, we take the trace over the degrees of 
freedom of the environment. The effective evolution, not 
necessarily unitary, is then given by: 



Hps)=Ti-e[Use(ps®\0) e (0\)uI e 
= Y, i1 e{iAU S e\Q) e Ps e 



U SE \P)E 



(6) 



where {\fj)} form an orthonormal basis for E, and the 
operator $ describes the evolution of the system S ( $ is 
usually called a quantum channel, in analogy with clas- 
sical communication theory [Tlj ) . Finally this evolution 
can be expressed only in terms of operators acting on S 
in the following form: 



$(ps)=J2 M nPsMl 



where the operators 

M„ = E (fj,\U S E\0) E 



(7) 



(8) 





U 




PS 





UpsTP 



a) 





Use 




Ps 




10) 





Hps 



b) 



FIG. 1: Unitary dynamics, a) Closed system b) Open system 
- $ describes the reduced evolution of S when we trace out the 
environment E. 



are the so-called Kraus operators [4y, [47|, [48| • The prop- 
erty J2p, M l M P = 11 guarantees that Tr[$(p s )] = 1, so 
that the operation $ is trace preserving. Furthermore, 
the evolution given by Eq. |7]) preserves the positive 
semi-definite character of ps - this means that $(ps) is 
also a density operator. It is important to note that the 
Kraus operators are not uniquely defined - performing 
the trace operation in Eq.© in different bases leads to 
different sets of equivalent operators, yielding different 
decompositions of the resulting density matrix. 

There are at most d 2 independent Kraus operators [n], 
[49l | . where d is the dimension of S. Together with 
Eq. (J8|), this property implies that, if {| is a basis 
in the space corresponding to S, then a dynamical evo- 
lution of S, corresponding to the Kraus operators {M^}, 
// = 0, . . . , cr — 1 , can be derived from a unitary evolution 
of S + E given by the following map: 



\<h) |o> 
\<h) |o> 



M !</»!) |0) 
M \cj> 2 ) |0) 



+ M d 2_ 1 U 2 ) 



J 2 
J 2 



\<j> d )\0) -> M |^)|0) + ---+Af d2 _ 1 | < /» d )|d 2 -l) , 

(9) 

where as before the operators M, act only on S. This 
map yields the guiding equations for our experiments. 

If the environment has many degrees of freedom (so 
that it can be considered a reservoir), then under Marko- 
vian and differentiability assumptions Eq. {7} yields a 
master equation (4c| . This is however less general than 
the Kraus approach, which applies even if the environ- 
ment has a small number of degrees of freedom. 



1. Global vs. Local environments 

If the system S is itself composed of N subsystems 
(Si , . . • , Sn), we must distinguish between two main 
types of environment: 

i) Global channels: in this case all the subsystems are 
embedded in the same environment, and can even com- 
municate through it. These channels perform non-local 
dynamics and, in principle, can increase the entangle- 
ment among the subsystems. 

ii) Local channels: each subsystem interacts with its own 
environment, no communication is present. The total 



evolution can be written as Us-lE-l 
Eq.© is replaced by: 



Us, 



and 



>M? PS M]? ®---®Mf . (10) 



This operation is clearly local, and therefore cannot in- 
crease the entanglement among the constituents. 

Obviously, for systems with N > 2 mixed dynamics is 
also possible, i.e, some subsystems interact with a com- 
mon environment and others with independent environ- 
ments. 
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FIG. 2: Monitoring the environment. The state of the environment 
is measured, post-selecting the state of S. 



2. Filtering operations - Monitoring the environment 

Instead of directly following the dynamics of system 
S 1 , one can infer it by monitoring its surroundings. For 
instance, by detecting a photon emitted by a two-level 
atom, we know for sure that the atom is in its ground 
state. This scheme is illustrated in Fig. O 

The formalism of the preceding sections must be 
changed to take into account the monitoring of the en- 
vironment. Rather than tracing over the environment, 
we perform a measurement on it. If the outcome i is 
obtained, the state of S evolves to: 



M iPs M\ 



Pi 



(11) 



where pi = Tr[MipsMj] is the probability of finding the 
outcome i. Notice that, if the state ps is initially pure, it 
will remain pure after the measurement on the environ- 
ment. 

This application of a single Kraus operator to the state 
is usually called a filtering operation [ll|j H3|- A se- 
quence of successive evolutions and measurements de- 
fines a quantum trajectory for the state of S - each 
record of the state of the environment defining a quantum 
jump [H[. 



III. COMPLEMENTARITY, QUANTUM 
CHANNELS, AND THE DYNAMICS OF 
ENTANGLEMENT 

Up to this point we were dealing with the open dynam- 
ics in a rather general way. From now on we specialize 
on systems composed of qubits, which are representative 
of many physical systems of interest for quantum infor- 
mation processing. Furthermore, we consider only local 
environments, which implies in the interaction of each 
qubit only with its own environment. This is the situa- 
tion for two decaying atoms separated by a distance much 
larger than the wavelength of the emitted radiation. The 
individual qubit dynamics is then used to describe how 
the initial entanglement of two qubits is degraded due to 
the action of these independent environments. 

An elegant way to illustrate how the entanglement of 
the system with the environment disturbs the individual 



properties of the subsystems is through the complemen- 
tarity relations presented in Ref. [3^|. This is described 
in the following sub-section. 



A. Complementarity relations 

A single qubit S in a pure state has two complementary 
aspects, particle-like, and wave-like [5l|, which is math- 
ematically expressed by the following relation [3?], H|| : 



Vj = 1, 



(12) 



where the ingredients are the single-particle predictabil- 
ity Vs and visibility Vs- The first is a measure of the 
single-qubit relative population, defined as = K ^)!- 
The second is a measure of single-qubit coherence and is 
defined as Vs = 2|(cr + )|. Here tjj, with i <G {x,y,z}, are 
the Pauli matrices, and a + = |1)(0|. 

When the qubit S gets entangled with an environment 
E, its state becomes mixed. This implies that another 
term should be included in the previous relation, which 
then turns into 13911: 



C 2 



SE 



(13) 



where Cse is the concurrence [18|, which measures the 
entanglement between S and E. Independently of the 
dimension of E, the bipartite concurrence of the pure 
composite state is defined as 



Cse = J2(l-Tr[p|]) 



(14) 



where p s = Ti e [pse] UM- 

We see from Eq. (| L3[) that whenever the entanglement 
between the two systems increases, the single-particle fea- 
tures are reduced. When Cse = 1, the visibility and pre- 
dictability vanish - the single-qubit state has then com- 
pletely decohered. 

Relation (|13p was tested experimentally, using nu- 
clear magnetic resonance techniques, in Ref. J52|. We 
present in section fVl experimental results for the dynam- 
ics of these three quantities, obtained with a linear optics 
setup. The complementarity relations among them will 
help us to understand the action of different types of en- 
vironments on qubits. 



B. Quantum channels 

We describe now some of the most usual channels for 
qubits; amplitude damping, dephasing, bit-flip, phase- 
flip, and bit-phasc-flip. 



1. Amplitude damping 

This channel represents the dissipative interaction be- 
tween the qubit and its environment. The emblematic 
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example is given by the spontaneous emission of a photon 
by a two-level atom into a zero-temperature environment 
of electromagnetic-held modes. 

A simple way to gain insight about this process is 
through the corresponding quantum map: 



\0)s\0)e 



\0) S \0) E ; 

V^pWsWe + Vp^sWe » 



(15) 



which can be traced back to the 1930 Weisskopf-Wigner 
treatment of spontaneous emission by an atom [53[ . 

The first line indicates that if no excitation is present 
in the system, it remains in the same state and the en- 
vironment is also untouched. The next line shows that 
when one excitation is present in the system, it can ei- 
ther remain there with probability (1 — p), or it can be 
transferred into the environment with a probability p. 

Notice that p in these equations is just a parameteri- 
zation of time. The relationship between the parameter 
p and time t for an atom interacting with an infinite 
number of electromagnetic field modes, initially in the 
vacuum state, under the Markov approximation, is given 
by p = (1 — e" rt ), where T is the decay rate. In this 
case, the state \1) E in the map above can be understood 
as one excitation distributed in all field modes. However, 
this map can also be used to describe the interaction of a 
two-level atom with a single mode of the electromagnetic 
field inside a high-quality cavity [54J. In this case the 
excitation oscillates between the atom and the field, and 
we should take p = sin 2 (fit/2), where fl is the vacuum 
Rabi frequency. 

The fact that the same set of equations describes the 
interaction with either a reservoir or an environment with 
a single degree of freedom is a consequence of the gen- 
eral character of the Kraus approach, as commented right 
after Eq. ©. 

These remarks show that it is actually very advanta- 
geous to describe the evolution of the system through a 
quantum channel, rather than through a specific master 
equation or Hamiltonian. Together with the parameter- 
ization of the evolution in terms of p, thus avoiding a 
specific time dependence, this leads to a very general de- 
scription, which includes many different processes in the 
same framework. 

In all cases, the dynamics represented by the map (|15[) 
has the following Kraus operators (in the computational 
basis {|0),|1)»: 



M = 







1 



VT - p 



VP 




(16) 



Let \x) = ol |0) + (3 |1) be a general initial qubit state, 
i.e, at p = 0. According to Eq. (0, it evolves under the 
amplitude channel to: 



H\X){X\) - I a *p y /TZT£ (i -p)\[3\ 2 



(17) 



We can see from this state that coherence decreases with 
increasing p. Also, the population of |1) is transferred 



to |0). When describing the spontaneous decay (p — 
1 — e~ r *), only in the asymptotic limit t — > oo coherence 
drops to zero and the system tends to the ground state. 
These conclusions can also be drawn from the expressions 
for the visibility and the predictability: 

T S (P) = ||a| 2 -|/3| 2 +2p|/?| 2 |Hl-2(l-p)|/?| 2 | ; 
V s (p) = 2y/T=p-\aP\ = VT - P Vs(0) , 

(18) 

where Vs(0) is the initial visibility. Furthermore, within 
the entire interval < p < 1 the qubit state is mixed. 
This is confirmed by the calculation of its entanglement 
with the environment: 



(19) 



Cs E (p) = 2\f3\ 2 Vp(l-p) 
which vanishes only at p = or p — 1. 

2. Dephasing 



Here the coherence of the qubit state disappears with- 
out any change in the populations. This process occurs 
often when a noisy field couples to a two- level system [55| . 

The corresponding unitary evolution map is given by: 



\0)s\°)e 
|1>s|0> b 



|0) s |0) £ , 

VT=P-\1) S \0) E + VP\1) S \1) E . [M) 



It can be understood as an clastic scattering, where the 
the state of the two-level system does not change, but 
the state of the environment undergoes a transition with- 
out any energy exchange, due for instance to the change 
of momentum of its constituent particles. Although the 
states of the computational basis {|0), |1)} do not change 
under this map, any superposition of them will get en- 
tangled with the environment. 

The characteristics of this type of channel can be ana- 
lyzed, as before, by observing the evolution of a general 
state. The corresponding Kraus operators are: 



M 



1 



Mi 







% /T - p J — v Vp . 

Therefore, the state \x) = a |0) + j3 |1) evolves to: 
H 2 af3*VT^P 



(21) 



(22) 



As previously stated, the populations do not change, 
as well as the state predictability: Vs{p) = | \a\ 2 — \(3\ 2 \ = 
■Ps(O). On the other hand, the visibility monotoni- 
cally decreases: Vs(p) = 2|«/?| >/l — p = \/l — p Vs(0), 
as the system S gets entangled with the environment 
R. The entanglement between them is easily evaluated: 
Cse(p) = 2\/p\ a f3\ = y/p^s(0). This emphasizes the fact 
that states with zero initial visibility do not get entangled 
with this type of environment. 
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TABLE I: Evolution of complementary aspects for the initial state \x) — ot |0) + j3 |1) under bit, phase, and bit-phase flip. 



Channel 




Vs{p) 




Vs(p) 


Csr(p) 


Bit flip 
Phase flip 
Bit-Phase flip 


(1 
(1 


-p)Vs(0) 
Ps(0) 
-p)Vs(P) 


1(2 
1(2 


- p)af3* + pa* f3\ 
(l-p)V s (0) 

— p)a(3* — p a* f3\ 


sjp{2-p)\a 2 -fi 2 \ 

Vp(2-p)Vs(0) 
^p(2-p)\a 2 + f3 2 \ 



3. Bit flip, phase flip, and bit-phase- flip 

In classical computation, the only error that can take 
place is the bit flip h 1. In quantum computation 
however, the possibility of superposition brings also the 
possibility of other errors besides the usual bit flip. They 
are the phase flip and the bit-phase flip. The first changes 
the phase of the state, and the latter combines phase- and 
bit-flip. 

The set of Kraus operators for each one of these chan- 
nels is given by: 

M = y/l-p/2 1 , Ml = 0i ; (23) 

where % = x give us the bit flip, i = z the phase flip, and 
i = y the phase-bit flip. These sets are easily interpreted 
as corresponding to a probability (1 — p/2) of remaining 
in the same state, and a probability p/2 of having an 
error. The factor of 2 in Eq. (f2"3")l guarantees that at 
p = 1 we have maximal ignorance about the occurrence of 
an error, and therefore minimum information about the 
state. The unitary maps for these channels are obtained 
by employing Eq. ((9]). 

In Table [H the evolution of the complementary as- 
pects, as previously defined, are summarized for these 
error channels. 



C. Entanglement dynamics 

Whenever the system S is composed of at least two 
subsystems, an initial entanglement among the subsys- 
tems evolves due to the interaction with the environ- 

mentfiE M, HE HE HE HE HE HE HE EuL HE HE EH 

I35l.l36|]. The detailed study of this process is of crucial im- 
portance for the implementation of quantum algorithms 
that rely on entanglement as a resource. 

Here we focus on two emblematic examples of entangle- 
ment evolution: the two qubit state \cf)) = a |00) + (3 |11) 
under local i) amplitude damping, and ii) dephasing 
channels. In the following analysis, the complementar- 
ity relation is not easily handled [52L HE, since it in- 
volves multipartite entanglement of mixed states. Nev- 
ertheless, in order to scrutinize the dynamics, we make 
use of similar figures of merit, namely: the bipartite vis- 
ibility (Vs 1 s 2 )i the concurrence between the subsystems 
(CsiS 2 )j an d the concurrence {Cse) between S = S\®S2 
and E = E\ (g) £/2- The definitions of these quantities 
follows. 



The bipartite visibility, 

V SlSa (p)=2K|ll)(00|)|, (24) 

measures the two-particle coherence for the state \4>) de- 
fined above. Notice that, given the initial state \(j>), and 
the fact that we are considering only local channels, this 
is the only coherence that plays a role in the dynamics. 

The initial pure state of the system S = Si ® S2 be- 
comes mixed when in contact with the environment. The 
degradation of the initial entanglement due to the cou- 
pling with the environment is quantified by the concur- 
rence defined in Ref. [l8l ]: 

C Sl s 2 (p)=max{0,A}, (25) 

where A = ^/\\ — s/M — V^3 — VM, with Ai's the eigen- 
values in decreasing order of: 

Ps 1 S 2 (p){v y <8> cr y )p* SlS2 (p)(cr y g) Oy) , (26) 

the conjugation being taken in the computational ba- 
sis {|00>, |01), 1 10), |11)}, and PSi s 2 (p) = $1 ® Hl&^l), 
where $i($2) is the channel applied to the first (second) 
qubit. 

The information spread from the initial pure state to 
the combined state - system plus environment - is related 
to the entanglement between S and E. The correspond- 
ing concurrence is fTsT ]: 

C SE (P) = ^(l-Tr^b)]) . (27) 

i) Amplitude damping - As described before in sec- 
tion llll BT1 the Kraus operators for this channel are given 
in Eq. (JTHJ). Under two identical local amplitude chan- 
nels, Eq. (fTUl) shows that the initial two-qubit state \(j>) 
evolves to the density operator 

/ \a\ 2 +p 2 \p\ 2 (l-p)a/T \ 

{1-p)p\I3\ 2 

(l-p)p\P\ 2 

V (l-p)a*/3 (1 -P) 2 |/3| 2 / 

(28) 

where the matrix is written in terms of the computational 
basis {|00), |01), |10), |11)}. 

The two-particle visibility is then Vs 1 s 2 (p) ~ 2(1 — 
p) I a/3 1 = (1 ^p)Vs 1 s 2 (0). The bipartite visibility decays 
linearly with p, reaching zero only when p = 1. 

For the entanglement between the subsystems, we 
have: 

C Sl s 2 (?) - max{0, 2(1 - p)\0\{\a\ - p\0\)} . (29) 
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p = 




FIG. 3: Two possible trajectories in the space of states un- 
der the action of amplitude damping, for initial states of the 
form a|00) + /3|11). The solid line represents a sudden-death 
trajectory, and the dashed line a case of infinite-time disen- 
tanglement. When p — 1, the two qubits are in the ground 
state. The border of the set of states is the locus of density 
matrices of incomplete rank. 



For the same initial concurrence (Cs 1 s 2 (ty = 2|ck/3|) , two 
entanglement decay regimes arc found: if \a\ > \(3\, then 
Cs 1 s 2 (p) > f° r all P G [0, 1), vanishing only at p — 1 
(as the visibility). However, for \a\ < \f3\ the entangle- 
ment between Si and S2 goes to zero at pesd = I <a / /3 1 
- the so-called entanglement sudden- death [32t [33j. If 
the parameterization (1 — p) = e~ rt is used, this implies 
a finite-time disentanglement, even though the bipartite 
coherence goes to zero only asymptotically. This phe- 
nomenon stresses that bipartite coherence is necessary 
for entanglement but does not coincide with it - the lat- 
ter being more fragile to noise. 

Entanglement sudden death requires the initial popu- 
lation of the doubly excited state |11) to be larger than 
the population of the unexcited state |00). This is re- 
lated to the fact that the state |11) is perturbed by the 
(zero temperature) environment, while the state 1 00) is 
not. Therefore, the bigger the initial "excited" compo- 
nent in \<j>), the stronger is the entanglement with the 
environment - thus leading to a faster decay of Cs 1 s 2 - 
Indeed, the entanglement between the system and the 
environment is given by: 



C se (p) = 2V2\p\y/p(l- P )y/l-\0Ml-p), (30) 

which increases when (3 increases and, for fixed /3, 
reaches its maximum for p = 1/2. This behavior is 
further stressed by realizing that the entanglement 
between each system and its own environment is also 
proportional to the excited-state amplitude: 



■-Si-E 



tip) = C S2 E 2 (p) 



Vp(i-p) ; (3i) 



vanishing only at p = and p = 1 . 

These two possible "trajectories" EH in the set of 
states are sketched in Fig. [3] For \a\ < \(3\ (solid line) the 
set of separable states is crossed at Pesd, thus the state 



becomes separable at finite time. However, for \a\ > \/3\ 
(dashed line), the state becomes separable only at p = 1, 
when the two qubits are in the ground state (|00) ) . 

This type of environment acts as a swapping process 
at p = 1, i.e., the state of the system (and the corre- 
sponding entanglement) is completely transferred to the 
environment 159, kjfj : 



(a|00)+/3|ll)) s ®|00) J 



0=1 



) s ®(a\OO)+0\U)) R . 

(32) 

The entanglement between the two environments is given 
by: 

C Ei e 2 (p) = max{0, 2p\0\ (\a\ - (l-p)\p\)}, (33) 

which shows that whenever there is entanglement sudden 
death for the two-qubit system, there is also sudden birth 
of entanglement (ESB) between the two corresponding 
environments (60| . The value of p for which ESB occurs, 
Pesb, is simply expressed in terms of the entanglement 
sudden death value: pesb = 1 — Pesd- This expression 
clearly shows that entanglement sudden birth may oc- 
cur before, simultaneously, or after entanglement sudden 
death, depending on whether pesd > 1/2, Pesd = 1/2, 
or Pesd < 1/2, respectively. 

ii) Dephasing - The evolved state in this case is given 
by: 



Pip) 



( H 2 (l-p)ap* 





V {l-p)a*P |/3| 2 



(34) 



As above, the two-particle visibility is given by 
V Sl s 2 ip) = 2(1 -p)\aP\ = (l-p)V SlSa (0), which leads 
to the same behavior as before. 

The entanglement between the subsystems is: 



C SlSa (p) = 2(l-p)|aj9| , 



(35) 



which is precisely equal to Vs x s 2 ip)- These two quantities 
have the same behavior as a function of p, and thus vanish 
at same point p = 1. There is no entanglement sudden- 
death. 

The entanglement between the system S and the envi- 
ronment E is given by: 



C SE ip) = 2\a(3Wpi2-p) 



(36) 



which reaches its maximum, for p fixed, when \a\ = |/?| = 
l/y/2. For every a and (3 the maximum of Cse as a func- 
tion of p is at p = 1, i.e, only when dephasing is com- 
pleted. However, increasing values of Cse do not imply 
sudden death of entanglement, since the corresponding 
state trajectory does not cross the region of separable 
states (see Fig. H}. 

One should note that, in contrast with the amplitude- 
damping case, here each system does not get entangled 
with its own environment: Cs 1 e 1 = C-s 2 E 2 = for 
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p = 




FIG. 4: Trajectory in the space of states under the action of 
dephasing, for initial states of the form a|00) + /3|11). The 
state is completely decohered at p = 1, when it reaches the 
borderline between entangled and separable states - and only 
then becomes separable. For this case, the trajectory stays 
always on the border of the set of states, since, for all p £ [0, 1] , 
the density matrix is not of complete rank. 



all p G [0,1]. This is expected, since the single-qubit 
visibility is zero at the beginning, and consequentially 
at all subsequent times (see Section IIIIB 2[) . What is 
though more surprising is that, apart from Cs 1 s 2 which 
monotonously decreases with p, all other two-quibt en- 
tanglements are identically zero. Decrease in the entan- 
glement of the two-qubit system is accompanied by the 
creation of legitimate multipartite entanglement. 
For p = 1, it is easy to see from Eq. (|2"0)) that 

H00) + /3|11)) S |00) B -> a|00) 5 |00) B + p\U) s \U) E , 

(37) 

which is a Greenberger-Hornc-Zcilingcr (GHZ) type [41( 
type of state, for which any two-qubit entanglement van- 
ishes. 

For arbitrary values of p, one can easily calculate 
the gen eralized multipartite concurrence proposed in 
Ref. gi| : 



C N = 2^ 2 /(2*-2)-£Tr(pf) 



(38) 



where the sum is over all nontrivial reduced density ma- 
trices of the N-particle system. We get 



Cs 1 s 2 e 1 e 2 (p) = '4 + 4p - p 2 
which monotonously increases with p. 



(39) 



IV. EXPERIMENTAL IMPLEMENTATION OF 
DECOHERENCE CHANNELS 

A. Single photons and multiple qubits 

Many experimental investigations of quantum infor- 
mation processes, such as basic quantum algorithms 



[6ll . |62| , |63[ , quantum teleportation [64( , and verification 
of new methods for measuring entanglement [65L [6t| have 
been based on the use of several degrees of freedom of sin- 
gle photons. While this type of a ppr oach does not lead 



to scalable quantum computation [67| , taking advantage 
of multiple degrees of freedom of photons allows for en- 
tan glem ent purification (r38j , improved Bell-state analysis 
f69ll7d.l7il .l72j and creation of high-dimensional entangle- 
ment [73j . The extra degrees of freedom have also been 
exploited to engineer mixed states through decoherence 

In the following we employ the polarization degree of 
freedom of a photon as the qubit, while its momentum 
degree of freedom is used as the environment. This choice 
enables us to implement controlled interactions between 
S and E. As in previous works [6l|, it is possible 
to implement a variety of operations on these two de- 
grees of freedom, using common optical elements such 
as wave plates and beam splitters. The formal corre- 
spondence between linear optics operations and one or 
multiple-qubit quantum operations has been provided in 
Ref. M- 



B. Sagnac Interferometer 



Fig. [5] a) shows a modified Sagnac interferometer that 
can be used to implement the dynamics discussed in sec- 
tion [Til An incident photon passes through a polarizing 
beam splitter (PBS), which splits the horizontal (H) and 
vertical (V) polarization components, causing them to 
propagate in opposite directions within the interferom- 
eter. The interferometer is aligned so that the 7J-path 
and y-path are slightly separated, which allows us to in- 
sert different optical elements separately into each path. 
The two paths then recombine at the same PBS, and 
are reflected or transmitted into modes or 1, depend- 
ing on the polarization. HWP(6># ) and HWP(<V) rotate 
the H and V polarization components of the incoming 
photon, respectively. If they are set at positions such 
that the polarizations are not rotated, the photon leaves 
the Sagnac interferometer in mode 0. If, however, a pho- 
ton, initially F-polarized, is rotated by HWP(#y) so that 
\V) — ► a \ V) + P\H), it will leave the interferometer in 
mode with probability |a| 2 , and in mode 1 with prob- 
ability |/3| 2 . The Sagnac arrangement is advantageous, 
since it is very robust against small mechanical fluctua- 
tions of the mirrors and polarizing beam splitter (photons 
in the two paths reflect off the same optical components) , 
as well as thermal fluctuations. The two optical paths are 
approximately the same, since they have identical lengths 
and both include a single half-wave plate. We now dis- 
cuss the implementation of decoherence channels with 
this interferometer. 
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a ) HWP(e H ) 




FIG. 5: Experimental apparatus for implementation of quan- 
tum maps and tomographic analysis, a) Sagnac interferome- 
ter. HWP(fly), HWP(fe) and HWP(0i) are half-wave plates 
aligned in such a way that 9y (8h, Oi) is measured with re- 
spect to V-polarization (H-polarization), PBS is a polarizing 
beam splitter, PP(</>) is a phase plate and unlabelled rectan- 
gles are mirrors. Polarization tomography is then performed 
on the output modes, b) Tomography setup used when trac- 
ing out the environment modes. The quarter- and half- wave 
plates QWP and HWP are used for the tomographic pro- 
cess and HWP (45°) is a half- wave plate used to incoherently 
recombine the environment modes at the PBS. c) Tomogra- 
phy setup used when monitoring the environment. Detectors 
DETo and DETi detect photons in modes and 1, respec- 
tively. 



1. Decoherence channels 

With the half- wave plates set to angles Oh and 9y, the 
Sagnac interferometer implements the transformation: 

\H) |0) — ► cos20 ff \H) |0) + sin20 H \V) |1) , (40a) 
\V) |0) — ► cos 29 v \V) |0) + sin20y \H) |1) . (40b) 

After the half-wave plate and the phase plate in output 
mode 1, the overall transformation is 

\H) |0) — >cos2d H \H) |0) +e^sin20 ff sin20i \H) |1) 
- e # sin26» i /cos26' 1 \V) |1) (41a) 

\V) |0) — >cqb29v\V) |0) +e # sin 29 v cos 26 x \H) |1) 
+ sin 20y sin 20i | V) |1) . (41b) 

By associating H and V polarizations respectively to 
the ground and excited states of the qubit, output modes 
and 1 to states of the environment, and adequately 
choosing the correct wave plate angles, a number of de- 
coherence channels can be implemented with this inter- 
ferometer. For example, setting 9h = 0, 9\ = 0, <j> — 
and identifying p = sin 2 20\/, the interferometer corre- 
sponds to the amplitude damping channel (|15p . 



TABLE II: Wave plate angles and phase 
herence channels. 


4> for different deco- 


Channel 


e H 


e v 


9i 


4> 


Amplitude decay 





e 








Phase decay 








tt/4 





Bit flip 


-6 











Phase flip 





-0 


tt/4 





Bit-Phase flip 


-6 


-0 





tt/2 



Using the same settings but with 0i = tt/4 implements 
the phase damping channel (f20|) . Also, the error channels 
shown in Table|T]can be implemented. For example, 9h = 
—9, 9y = 9, 9\ = <j) = implements a bit-flip channel 
with p = sin 29 v . Table HJ shows the wave plate settings 
for several different decoherence channels. 

In order to investigate decoherence using this interfer- 
ometer, it is necessary to combine modes and 1 incoher- 
ently for detection, which is the experimental equivalent 
to mathematically "tracing out" the environment. This 
was done using the two input ports of the PBS used for 
polarization tomography, as shown in Fig. [5] b). In- 
coherent combination is guaranteed as long as the path 
length difference is larger than the coherence length of 
the photon. Using two or more interferometers, one can 
use similar setups to study the evolution of multipartite 
states subject to different combinations of independent 
error channels. 



2. Monitoring the environment 

As opposed to uncontrolled physical processes that in- 
duce decoherence in other systems, our interferometric 
arrangement allows us to monitor the environment. This 
can be done for instance by simply detecting photons 
in momentum modes or 1 individually, as illustrated 
in Fig. [5]c). With this setup we are able to experimen- 
tally investigate the filtering operations and the quantum 
jumps described in section III A 21 The corresponding ex- 
perimental results are discussed in the next section. 



V. EXPERIMENTAL RESULTS: SINGLE-QUBIT 
DECAY AND THE DYNAMICS OF 
COMPLEMENTARITY RELATIONS 

In the experiments reported in this section and the 
next, we controlled the system-environment interaction 
by varying the parameter p. For each value, we per- 
formed full tomography of the single or two-photon polar- 
ization state and reconstructed the density matrix using 
the maximum likelihood method [l^, [78| . The purity and 
the concurrence were obtained from the reconstructed 
density matrix unless otherwise noted. The theoretical 
predictions were obtained by evolving the reconstructed 
initial state, corresponding to p = 0, using the Kraus 
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operator formalism discussed in section [TTJ Vertical ex- 
perimental error bars were determined by Monte-Carlo 
simulation of experimental runs obeying the same Pois- 
sonian count statistics. The value of p was determined 
by one of two methods. In our earlier experiments, we 
used the direct readout of the angle (larger horizontal 
error bars, due to coarse angular setting). In later exper- 
iments, we improved on this by developing a simple way 
to determine p empirically, which we now describe. First, 
we block the interferometer arm corresponding to the H 
component [propagation in the counterclockwise direc- 
tion in Fig. [5] a)], and measure the counts Co in output 
mode 0, with the tomographic plates set for measuring 
V polarization. Then, still blocking the H interferom- 
eter arm, we measure the counts c\ in output mode 1, 
with the tomography system set to H. We obtain p from 
p = ci/(co + c{). This method is more precise, since the 
uncertainty in p comes from photon count statistics. 



A. Amplitude damping channel 




0.2 0.4 0.6 0.8 1.0 



FIG. 6: (Color online) Evolution of the quantities involved in 
the qubit-environment complementarity relation (|13[) under 
the amplitude damping channel: V% (blue circles), Vf (red 
squares) and C% E (green diamonds). Also shown is the sum 
V% + Vf + Cg E (yellow triangles). The solid lines are the 
corresponding theoretical predictions, obtained by applying 
the amplitude damping channel to the experimentally recon- 
structed initial state. 



For the study of the amplitude decay of a single qubit, 
we use a c.w. solid state laser (405 nm) to pump a 5 mm 
long LUO3 non linear crystal, producing photon pairs 
from spontaneous parametric down conversion (SPDC). 
The signal and idler photons are prepared in polariza- 
tion product states with V polarization. Here the idler 
photon is used only as a trigger, and is sent directly to a 
detector equipped with an interference filter centered at 
800nm, (65nm FWHM) and a 0.5mm diameter pinhole. 
The signal photon goes through the interferometer shown 
in Fig. [5] a) , with wave plates aligned for implementation 
of the amplitude damping channel, as discussed in section 
IIVB 11 After the interferometer, modes and 1 propa- 
gate through wave plates and a polarizing beam split- 
ter, necessary in the tomography process. Afterwards 
they are detected through an interference filter centered 
around 800nm, with lOnm bandwidth and a 0.5mm pin- 
hole. Coincidence counts arc registered with counting 
electronics and a computer. 

The amplitude damping channel was implemented for 
a single qubit, with the detection system set to trace over 
the modes of the environment, using the detection setup 
shown in Fig. [5]b). The input polarization was prepared 
in a superposition state a \H) + (3 \V), with \/3\ > \a\. 

It is illustrative to view the effect of the channel by 
measuring the quantities involved in the complementarity 
relation discussed in section IIII Al In Fig. [6] we show the 
evolution of the squared predictability Vg, the squared 
visibility V|, and the system-environment entanglement, 
quantified by the squared concurrence, C SE , as a function 
of p for the same initial state as above. The concurrence 
Cse was calculated from the density matrix using (|14[) . 
and coincides with what is expected from Eq. (fT9"l) . Vs 
and Vs were determined directly from the polarization 



measurements using 



c H + c v 

V S = J(^_- l)' ,43) 
y \CH + c v J \CH + cv J 

for each value of p, where Cj are the number of counts 
with j polarization, with + and R corresponding to 45° 
linear polarization and right circular polarization, respec- 
tively. It can be seen from Fig. [6] that both quantities 
agree with Eqs. lfT5]). Though V%, Vf , and C 2 SE evolve 
with p, the sum of these three quantities (yellow triangles 
in Fig. [5]) satisfies relation (flU)) for all p. 



B. Monitoring the environment 

We demonstrate now a peculiar effect of the dynamics 
of open quantum systems. If the qubit, under the ac- 
tion of the amplitude-damping channel, is initially in a 
superposition of the states |0) and |1), and we monitor 
the state of the environment, finding it with no excita- 
tion at all times, we still observe a decay of the system 
towards the ground state. This can be understood as fol- 
lows: even if there is no energy transfer between system 
and environment, by constantly monitoring the environ- 
ment and finding no excitations in it, we gain information 
about the system, which is expressed as a change in its 
state. For example, consider the arrangement used to 
implement the amplitude damping channel (|15[) . with an 
initial state (a \H) + 0\V))s <8> |0) £ . This state evolves 
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P 

FIG. 7: (Color online) Evolution of the population V when 
monitoring (red squares) or tracing out (blue circles) the en- 
vironment. Error bars are unnoticeable in this scale. The 
lines are the corresponding theoretical predictions (red solid 
line and blue dashed line). 
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FIG. 8: (Color online) Evolution of the purity of the qubit 
state when monitoring (red squares) or tracing out (blue cir- 
cles) the environment. Error bars are unnoticeable in this 
scale. The lines are the corresponding theoretical predictions 
(red solid line and blue dashed line). 



to 



\^(p))=a\H)\0)+(3^T^\V)\0)+/3^\H)\l). (44) 
Tracing out the environment, the polarization state is 



q(p) 



\a\ 2 + \(3\ 2 p af3*^T~^ 
a*(3VT=p- |/?| 2 (1-P) 



(45) 



with p 2 = sin2#y, whereas, projecting onto the "unex- 
cited" |0) state of the environment, the polarization state 
becomes 



a\H)+f3^T^p-\V) 

N 2 + |/3| 2 (i-p)] 1/2; 



(46) 



which "decays" to \tp(p = 1)) = \H), just as g(p = 1) 
given in (j45]l . We illustrate this phenomenon by com- 
paring the dynamics of a qubit under the action of the 
amplitude damping map (|15[) . for two cases: (i) when 
we trace out the environment's degrees of freedom (using 
the tomographic setup illustrated in Fig. E]b), and (it) 
when we monitor its state in the unexcited state, using 
the tomographic setup shown in Fig.[5]c). As above, the 
input polarization is prepared in a superposition state 
a\H)+p\V), with |/3| > \a\. 

Figure [7] shows the evolution of the population V, for 
both cases. We see that not only are the two dynamics 
different, but also that the decay takes place even if no 
excitations are transferred to the environment. When the 
environment is traced out, the linear evolution is equiv- 
alent to an exponential decay (linear in p), while in the 
case where the environment is monitored the decay is 
retarded. 

Figure [5] shows the evolution of the purity in these two 
cases. We see that when we monitor the environment, 
the system is always close to a pure state. The little 
mixedness arises from the fact that our initial state is 
not perfectly pure. 



In these figures, the lines are the theoretical predic- 
tions, which are obtained (as introduced in section lTl A2p 
by using the Kraus operators corresponding to the am- 
plitude damping channel (fT5")) . When tracing out the 
environment, both operators Mq and Mi are used, while 
when monitoring the environment in its unexcited state, 
only the no-jump operator Mq is used, the resulting state 
being renormalized afterwards. The agreement between 
theory and experimental data is quite good. These re- 
sults show that, even though no excitation is transferred 
to the environment, the continuous acquisition of infor- 
mation about this fact changes the state of the qubit, 
increasing the probability that it is found in the unex- 
cited state. 

This phenomenon allows the distillation of entangle- 
ment of a two-qubit system through continuous local 
monitoring of the corresponding independent environ- 
ments. Indeed, for an initial state a|00) + /3|H), with 
\a\ < |/3|, continuous monitoring of the unexcited en- 
vironment leads to increase of the 1 00) component, im- 
plying that the state approaches a maximally entangled 
state, before decaying to the state |00). Within the 
framework of the Sagnac- interferometer setup, applied 
to each of the two entangled photons (as described in 
detail in the next Section), the evolution of the system 
under continuous monitoring corresponds to measuring 
both qubits in output mode 0, while p changes from 
to 1. If the two-qubit state is |*(0)) = a\HH) + /3\VV) 
for p = 0, then the state for p ^ 0, conditioned to the 
measurement of output mode 1 in |0) for both qubits, is 
given by: 



I*(P)>C 



a\HH) + (3{l-p) \VV) 

[|«| 2 + |/3| 2 (1 -p) 2 ]l/2 



(47) 



Setting 1 — p = \a/j3\ yields a maximally entangled state. 

Therefore, continuous monitoring of the environment 
for the two-qubit case corresponds to a quantum distil- 
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FIG. 9: (Color online) Experimental setup for experimental 
investigation of entanglement dynamics under decoherence. 
Polarization-entangled photons pairs were created using the 
two-crystal SPDC source [§(}, pumped by a 442 nm c. w. He- 
Cd laser. With proper spatial and spectral filtering, the 884 
nm photons created in the two nonlinear crystals are prepared 
in an entangled state of the form |9> = \a\\HH) + \f3\e iS \VV). 
The coefficients \a\ and \j3\, and the relative phase S are 
controlled by manipulation of the pump-beam polarization, 
which is easily realized with a half-wave plate (HWP) and a 
tilted quarter- wave plate (QWP) [§(}. Tomographic analysis 
was performed with additional HWP's, QWP's and polarizing 
beam splitters (PBS). 



lation scheme [7£ 







FIG. 10: (Color online) Tomographic reconstruction of the 
real part of the density matrix for different values of p for an 
initial state close to |01) given in Eq. (|48ap . 



its evolution through the relevant Kraus operators. The 
dynamics of entanglement was investigated under the ef- 
fect of two different decoherence channels, implemented 
by sending the twin photons through independent Sagnac 
interferometers. Full bipartite polarization-state tomog- 
raphy [77], [78[ was performed for different values of p. 



VI. EXPERIMENTAL RESULTS: THE 
DYNAMICS OF ENTANGLEMENT 

Using two Sagnac interferometers, we studied the dy- 
namical behavior of global properties of an entangled pair 
of photons generated with spontaneous parametric down 
conversion. Fig. [9] shows the experimental setup. 

The source was arranged to generate pairs of photons 
in one of two non-maximally entangled states given by 

[01> = \\HH) + ^e ie \VV), (48a) 
\®2) = ^-\HH) + \e iS \VV). (48b) 

These states contain the same amount of entangle- 
ment: the inital concurrence is ideally C = 2\a0\ = 
V3/2 ~ 0.87. However, we measured C = 0.82±0.04 and 
C = 0.79 ± 0.11, respectively, due to the fact that they 
were not 100% pure. The decrease in purity is mostly 
due to small imperfections in the mode matching of the 
interferometers, and the ang ular dependence of the phase 
of the two-photon state j80[. To simplify the description, 
we refer to the initial state in the following as either (|48a|) 
or (|48b[) . meaning in fact that it was close to these states. 
The theoretical predictions are derived from the actual 
experimental state obtained when p = 0, by calculating 



A. Amplitude damping - Entanglement sudden 
death and entanglement witness 

Using the wave plate configurations listed in Table [Til 
dual amplitude damping maps (|15[) with the same p were 
implemented for the initial state (|48ap . Tomographic re- 
constructions of the real part of the density matrix for 
different values of p are shown in Fig. [TDJ The cor- 
responding analytical expression is given by Eq. (|28[) . 
They illustrate the evolution of the populations and co- 
herences as a function of the parameter p. 

Figure[TT]displays the experimental results for the con- 
currence (|2"5)) . The theoretical prediction (denoted by the 
full line in the figure) was obtained by applying Eq. (|25|) 
to the evolved density matrix, which in turn is deter- 
mined by applying the Kraus operators (|16|) to the re- 
constructed density matrix for p = 0. The vanishing 
of the entanglement for a transition probability p < 1, 
corresponding to a finite time, is clearly demonstrated 
in Fig. 111! This phenomenon was termed entanglement 
sudden death [32| . and our setup allowed for its first ob- 
servation, which was previously reported in [33| . 

Also shown in Fig. |TT| are the results obtained from an 
entanglement witness, evaluated at each data point. An 
operator W is an entanglement witness if Tr(Wp) > 
for any separable state, and there exist entangled states 
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FIG. 11: (Color online) Global properties under the ampli- 
tude damping channel for the state (|48a[) : A (as defined af- 
ter Eq. (|25p . red circles) and the entanglement witness (|52p 
(green squares) as functions of the transition probability p. 
The theoretical prediction for the concurrence (red solid line) 
is derived from the density matrix obtained by applying the 
amplitude decay channel to the experimentally reconstructed 
initial state. Disentanglement occurs for p < 1. 
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FIG. 12: (Color online) Global properties under the ampli- 
tude damping channel for the state (|48b|) : A (as defined af- 
ter Eq. (|25[) . red circles) and the entanglement witness (|52p 
(green squares) as functions of the transition probability p. 
The theoretical prediction for the concurrence (red solid line) 
is derived from the density matrix obtained by applying the 
amplitude decay channel to the experimentally reconstructed 
initial state. For this state, disentanglement occurs only at 
p = l. 



a for which Tr(Wa) < 0. For initial states of the form 
\a\\HH) + \0\exp(iOo)\VV) (49) 

and the amplitude decay channel, it is possible to define 
a "perfect" p- independent witness [3l|, so that — Ti (Wp) 
coincides precisely with A in Eq. ([2"5]) , thus yielding the 
concurrence for all p. It is given by 



W flo = l-2|#(0 O )> ($(0o)| , 



where 



1*0?)) 



V2 



(\HH) 



Then it is easy to show that 



-Tr 



We p(t) 



\VV)), 



(50) 



(51) 



(52) 



The concurrence 



where P(9,t) = Tr[|$(0))($(0)|p(i)]. 
?5l) can be written as 

C[p(t)) =max{0,r eo } 



That is, the concurrence is equal to twice the excess prob- 
ability (with respect to 1/2) of projecting the system in 
the maximally entangled state |$ (0o))- It is remarkable 
that in this case the concurrence can be given a simple 
physical interpretation, and moreover that this is valid 
throughout the evolution of the system (which means, in 
our case, that it is independent of p). The concurrence 
could then be determined directly by measuring the prob- 
ability of finding the system in the maximally- entangled 
state |$ (0 O )). 



In our experiment, however, the initial state is not 
pure, so We is not a perfect witness. In order to com- 
pute the best witness in this case (which yields the up- 
per bound of — Tr[Wep]), we choose theta for each data 
point as the argument of the pvvhh element of the re- 
constructed density matrix, and then obtain T through 
Eq. ([52")) . The same could be achieved by projecting the 
state of the system on | x &(0)) and scanning in order to 
get the minimum value for — Tx[We/9(t)]. 

One should note that the initial phase 0o is not changed 
by the amplitude damping channel, as shown by Eq. (|28|). 
so in principle this procedure should be adopted only for 
the p = state, the corresponding witness being then 
valid for all values of < p < 1. However, in the 
experiment, changing the angle of the half- wave plate 
HWP(0y) in Fig. O actually affects the corresponding 
optical path, due to imperfect alignment of this plate, so 
a p-dependent phase shows up between the states | V) |0) 
and |if)|0) on the right-hand side of Eq. |40b|) . This 
phase does not affect the concurrence, but it implies that 
the best witness depends on p. For this reason we find 
for each data point, from the reconstructed density ma- 
trix. Fig. [11] shows that in this case T underestimates 
the entanglement. 

For the state |62) denned in Eq. (|48b|) the situation 
is drastically different. Fig. [12] shows the concurrence as 
a function of p. In this case the entanglement disappears 
only when p = 1. Also shown is the entanglement witness 
([52]) calculated from the reconstructed density matrices, 
which is always less than the actual value of the concur- 
rence. The witness is not optimal since the initial state 
is not completely pure. 

Figures [TT1 and [T2l together constitute an experimental 



14 



1.0 



0.8 




0.0 0.2 0.4 0.6 0.8 1.0 



P 

FIG. 13: (Color online) Concurrence for the initial state (|48a[) 
under the phase-damping channel: the circles correspond to 
the concurrence determined from the reconstructed density 
matrix for each value of p, while the line is the correspond- 
ing theoretical prediction, obtained by applying the phase- 
damping channel to the initial (p — 0) density matrix. Dis- 
entanglement occurs only at p = 1. 
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FIG. 14: (Color online) Evolution of the purity for the dif- 
ferent decoherence models: amplitude damping (red circles) 
and phase damping (blue squares), and their corresponding 
theoretical predictions. 



VII. CONCLUSIONS 



confirmation that two states with the same initial amount 
of entanglement may follow different decoherence "trajec- 
tories" in the space of states, as discussed in section lHICI 
and illustrated in Fig. [3] 



B. Phase damping 

As shown in Table [Til by adjusting the wave plate 
angles the Sagnac interferometers implement the phase 
damping channel (|20[) . Experimental results for the con- 
currence arc presented in Fig. [13] for the initial state 
(|48a| . There is no sudden death of entanglement, and 
concurrence vanishes only when p = 1. 



C. Evolution of purity 

Figure[14]shows the evolution of the purity of the initial 
state (|48a|) for the amplitude damping channel (fT5|) and 
the phase damping channel |20|) . For the phase damping 
channel, the change in the purity is monotonous. For 
amplitude damping there is an initial decrease, and then 
it increases again up to 1 for p = 1. The difference in 
behavior reflects the fact that amplitude damping pro- 
motes a swapping between system and environment, so 
that the system ends up in state \HH), while dephasing 
leads to an increase of multipartite entanglement, with 
the system plus environment evolving towards a GHZ- 
likc state. 



We have presented a scries of experiments that inves- 
tigate the dynamics of entangled open quantum systems, 
and also the dynamics of a single qubit under continu- 
ous monitoring of the environment. By adjusting a set 
of wave-plates, our linear-optics setup is capable of im- 
plementing a number of singlc-qubit decoherence chan- 
nels. We present experimental results for the amplitude- 
damping and phase-damping channels for single and two- 
qubit systems. Decoherence of a single qubit is investi- 
gated through the use of complementarity relations. The 
effect of decoherence on entanglement , including the phe- 
nomenon of entanglement sudden death, is experimen- 
tally demonstrated. 

Our setup has an appealing feature: it allows the inves- 
tigation of filtering operations, implemented by monitor- 
ing the environment. This is an experimental realization 
of quantum trajectories [44| . which lead to a description 
of the interaction of a system with an environment in 
terms of pure states. For amplitude damping, our ex- 
perimental results demonstrate that it is possible to in- 
duce decay of a system by verifying, through continuous 
measurements, that no excitation is transferred to the 
environment. We have shown that this procedure, for an 
initially two-qubit entangled state, is equivalent to en- 
tanglement distillation. 

The experimental investigation of the environment- 
induced decay of entanglement in other systems (see, for 
instance, Ref. (35j ) is of course of fundamental impor- 
tance, and should help to throw new light on the subtle 
relation between local and global dynamics of entangled 
systems. 

The parametrization of the quantum channels consid- 
ered in this paper, in terms of the transition probability 
rather than time, is very convenient. It accommodates 
different kinds of dynamical behavior in a universal de- 
scription, which allows one to extend the realm of appli- 
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cation of the results obtained here: they include not only 
the decay of two-level systems interacting with individual 
and independent environments, but also the oscillatory 
exchange of energy between each qubit and another two- 
level system, which could be for instance the vacuum and 
one-photon subspace of a cavity mode. 

This is a quite advantageous strategy for investigating 
the dynamics of disentanglement, since it pinpoints the 
main features of this process within a quite encompassing 
framework. In fact, rather than the investigation of a 
particular physical system, our procedure amounts to the 
experimental implementation of quantum maps which, 
due to their generality, play a very fundamental role in 
quantum information. 
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